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ABSTRACT 

We present high resolution NICMOS images of random fields obtained in 
parallel to other HST observations. We present galaxy number counts reaching 
iJi.e = 24. The ifi.e—band galaxy counts show good agreement with the deepest 
I— and A—band counts obtained from ground-based data. We present the 
distribution of galaxies with morphological type to ifi.g < 23. We hnd relatively 
fewer irregular galaxies compared to an /—band sample from the Hubble Deep 
Field, which we attribute to their blue color, rather than to morphological 
K-corrections. We conclude that the irregulars are intrinsically faint blue 
galaxies at z<l. 


Subject headings: cosmology: observations — galaxies: evolution — galaxies: 
fundamental parameters — galaxies: irregular — infrared: galaxies — surveys 


1. Introduction 

The Near-Infrared Camera and Multiobject Spectrometer (NICMOS; Thompson et 
ah 1998) has taken multi-color images and grism spectra of random fields in the NICMOS 
Parallel Program (c/. Holfeltz et al. 1998). Like the Medium Deep Survey (MDS; Griffiths 
et al. 1994) and the Hubble Deep Field (HDF; Williams et al. 1996), the NICMOS parallel 
imaging provides high-resolution images of galaxies at the faintest magnitudes observed. 
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Our understanding of the evolution of galaxies has largely depended on our knowledge 
of their place in the Hubble sequence. The MDS provided the hrst measure of faint 
galaxy number-magnitude counts as a function of morphological type (Driver et ah 1995, 
Glazebrook et ah 1995). It showed a population of rapidly evolving irregular/peculiar 
galaxies that appear to explain the faint blue galaxy excess. At large redshifts, it becomes 
problematic to deduce a galaxy’s morphology from optical images, as the observation 
becomes increasingly dominated by blue and UV light from young stars, leading to a 
“morphological K-correction” (O’Connell 1997). With NICMOS deep imaging, it is now 
possible to sample the older stellar population and measure unbiased galaxy morphologies. 

NICMOS parallel imaging also provides galaxy number counts, as deep or deeper than 
are possible from the ground. Both optical and infrared number count surveys have detected 
excesses of faint galaxies over simple no-evolution model predictions (Tyson 1988, Gardner 
et al. 1993; see Ellis 1997 for a review). Theoretical models can explain the continued rise 
of the K-band counts, seen in the recent surveys from the Keck telescope, using either 
evolutionary effects or a dwarf-rich luminosity function - that is, one with a steep faint-end 
slope (c.f. Bershady et al. 1998). However, ground-based counts at the K=23 level are 
based on only ~ 50' total over 6-7 helds of view, so the critical last few magnitude bins 
may be affected by biased lines of sight. 

In this paper we present an analysis of Hi^ images obtained in the NICMOS parallel 
program in 1997. In section 2 we describe the data and our reduction methods. In section 3 
we present the galaxy number counts; in section 4 we discuss a quantitative morphological 
classihcation of the galaxies and compare the result to optical surveys. 


2. Archival Data and Data Reduction 

The parallel program consisted of “JHK” imaging using cameras 1 and 2 (hereafter 
NICl and NIC2) at their best focus until November 1997, and Ji.i and Hi^ imaging in 
Camera 3 (NIC3), as well as some grism spectroscopy, thereafter. Ri.g data were taken in 
individual exposures of 256s or 1024s for NIC2 and 256 or 576 seconds for NIC3. We selected 
160 NIC2 pointings at galactic latitudes above 20° that are unaffected by foreground stars 
or the primary target and have no instrumental oddities or other impediments (a total 
area of 15.2 square arcminutes). We also examined the hrst 30 NIC3 helds (20.2 square 
arcminutes). An independent analysis of the NIC3 data is presented by Yan et al. (1998). 

NICMOS parallel data are taken in “MultiAccum” mode, in which the array is read 
out nondestructively many times during each exposure. The resulting data cube contains 
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multiple measurements of the accumulated flux at increasing integration times. The count 
rate for each pixel can be determined from the slope of a line fit to these data (after dark 
subtraction). This procedure allows cosmic ray rejection by identifying breaks in the slope. 

NIC3 data were reduced using the CALNICA IRAF routines (Bushouse, 1997), using 
a model (noiseless) cube of multiaccum dark data constructed by STScI and a skyflat 
derived from the median of all NIC3 Hi q observations available at the time. NIC2 data 
were reduced independently of the CALNIC routines which were still under development at 
the time. The major difference between our routines and CALNICA is that we £t the slope 
of the multiaccum data before flatfielding. There was insufficient sky background in the 
NIC2 data to construct a viable skyflat and the best inflight pointed flat was used instead. 

Individual reduced frames were registered for each held. For NIC3 data, we used 
the centroid of many objects to calculate the precise offsets between dithered frames. 
Registration was performed with bilinear interpolation for fractional pixel shifts. During 
registration, a second cosmic ray rejection was performed using a sigma-clipping algorithm. 
Many NIC2 fields contained no bright objects for registration. In that case, we used the 
astrometry information from the telescope pointing to calculate fractional pixel shifts. 

Objects were identified using the SExtractor software (Bertin & Arnouts 1996), 
admitting no object with fewer than 6 pixels. Identihcations were carefully checked by 
eye to avoid false detections from cosmic rays or noise at the vignetted edges of NIC3 
frames. Photometry was performed in elliptical apertures set at 2.5 times the Kron “first 
moment” radius (Kron 1980). We measured Icr limiting magnitudes by taking the standard 
deviation of aperture photometry performed on several hundred randomly selected positions 
throughout the frame. In a typical photometric aperture (3.5 pixel radius) we measured 
galaxies with Hi 6=22.8 for NIC2 (and 22.9 for NIC3) at the 3a level in 1024 second 
exposures. Limits are expressed in magnitudes on the Vega system. The offset to the 
AB magnitude system (Oke 1974) is -1-1.3 for the iJi.e Liter. The deepest Lelds reach 
Hi 6 — 24.3 in NIC2 and 24.0 in NIC3. In the 160 NIC2 Lelds, 142 extended objects were 
detected. In the 30 NIC3 Lelds examined, more than 500 objects were detected. 


3. Galaxy Counts 

We examine the galaxy content of the central “overlap” regions of the registered 
images, that is the regions with full exposure time. We estimate completeness with 
simulations based on recovery of mock objects inserted into real coadded images. We £nd 
> 90% completeness at the 5a level and 60% completeness at 3a. We have included a small 
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(~ 10%) correction to our completeness estimate to account for reduced sensitivity to larger 
objects (see Bershady et al. 1998). For objects bright enough for a reliable estimate of 
the full-width half-maximum (based on the SExtractor algorithm), we separate unresolved 
stars from the galaxies. In a typical held, this limit is Hi.e ~ 22.5. We then assume that all 
fainter objects are galaxies, since there should be little contamination from faint stars in 
these helds. A similar assumption was made by Djorgovski et al. (1995) in making the deep 
K-band galaxy counts. 

Figure 1 compares the band counts to the deepest ground-based K-band counts, 
using a model prediction of their median color (the ncmod program, Gardner 1998) to 
translate the counts from the K-£lter to the NICMOS Hi^ filter. For example, ncmod 
predicts Hi q—K = 0.91 at Hi q=23. As has been seen in K-band data (Bershady et 
al. 1998, Moustakas et al. 1997), the Tfi.e galaxy counts continue to rise at the faintest 
levels. We note that the NIC3 counts appear less smooth than individual ground-based 
surveys or the NIC2 counts. While this effect could be real, it may be the result of 
non-Poissonian errors (e.g. held to held variations, completeness,etc.) Our future analysis 
of the entire NIC3 sample will fully explore this question. Our results show that NICMOS 
has already reached Hi.g ~ 24.3 at the 5cr completeness limit, using the most conservative 
error estimates. In the longest future pointings, NICMOS parallel data will go deeper than 
current ground-based IR number counts, and will accumulate more area. 


4. Morphology 

We take advantage of the high resolution of NICMOS images to evaluate galaxy 
morphology using the classification algorithm of Abraham et al. (1996a). The procedure 
uses the central concentration (C) and asymmetry (A), to separate E/SO’s, spirals, and 
irregulars. Ellipticals have strong central concentration, and low asymmetry. Spirals show a 
considerable range of asymmetry and lower central concentration. Irregulars, peculiars and 
mergers he in the “exclusion zone” of asymmetric and diffuse objects, where no ellipticals 
or spirals are expected to fall. 

The asymmetry is determined by rotating the galaxy image about its center, subtracting 
the original and measuring the ratio of hux in the absolute value of the difference to that 
in the original. The central concentration compares the hux in an aperture containing the 
inner third of the object as dehned by the second order moments to the total hux (Abraham 
et al. 1994). Both C and A are very sensitive to the size and center of the aperture 
containing the object. Abraham et al. dehned the center of the object to be the brightest 
pixel in the object after convolution with a 1 pixel wide gaussian. They dehned the object 
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aperture to be an isophot enclosed by the pixels l.Scr above the sky. We have used the same 
algorithm to determine C and A, but we use the object centroid as the center, and we use 
an elliptical aperture determined by SExtractor, also based on l.Scr sensitivity. 

To calibrate our classihcation routines, we have used a reference sample of 81 
ground-based images of E through Scd galaxies obtained by Frei et ah (1996) from Lowell 
observatory. This is the same comparison set used by Abraham et ah Comparing the 
values of C and A we measure for this sample to the Abraham et ah values, we obtain 
similar results with a scatter roughly equal to the errors, cr^ ~ crc ~ ±0.05. We also obtain 
consistent C and A values for artificially redshifted versions of the Frei sample. Abraham et 
ah applied a complex artificial redshifting scheme, applying a K-correction pixel by pixel. 
We have instead concentrated on understanding the effects of the pixel scale. We hnd that, 
as expected, galaxies can appear more symmetric at lower resolution, while their central 
concentration is almost unchanged, thereby preventing a spurious increase in the number 
of irregular galaxies when classihed at NIC2 resolution. To conhrm this expectation, we 
artihcially redshifted the Frei sample to z = 0.3, 0.5, and 0.7 with the NICMOS plate 
scale. We continue to see a clean separation of morphological types agreeing with the 
unredshifted type. At higher z, some galaxies cross the Elliptical/Spiral boundary in about 
equal measure in each direction, but we see no suggestion that changing the sampling or 
signal-to-noise of galaxies moves them into the “irregular/peculiar” part of the plane. 

We automatically classihed all galaxies with ^lOcr, the cutoff at which our “eyeball” 
classihcations could be trusted (Hi.e < 23 in the deepest helds and 1 to 2 magnitudes 
brighter in the shallowest). Figure 2 shows the morphological classihcation of all such 
objects in our NIC2 number counts sample (106 galaxies total). We hnd that based on our 
visual identihcation, the Spiral/Elliptical diherence could be identihed mostly using central 
concentration, while Spiral/Irregular separation relies heavily on asymmetry. However, no 
fewer “misidentihcations” would be achieved using a vertical cut in the C/A plane, so we 
retain the divisions of Abraham et ah, as borne out in our simulations. 

We hnd few irregular galaxies (17% of the total sample and 23% of galaxies fainter than 
77i.6=21). By comparison, the MDS survey found that 21% of the galaxies at /814 < 22 are 
irregular, while in the HDF the percentage rises to 30-40% at 22< /814 < 25 (Abraham et 
ah 1996b). Our i7i.6=20-23 limits correspond to 1814 = 22 . 1 - 24.5 for a mixture of elliptical 
and spiral galaxies, or /8i4=21.3-24 for late type spirals and irregulars. This quantitative 
result is conhrmed by visual inspection; no “train wrecks” are seen (see Figure 3). 

This fraction of faint, infrared-selected irregular galaxies appears small compared to 
the number found in the MDS and HDF. There are two possible explanations for this ehect: 
1) We are seeing the same galaxies but they are more regular in the NIR, or 2) we do not 
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see all of the irregulars because they are too blue. Put another way, are we seeing different 
stars or different galaxies? 

To distinguish between these cases, we compare our results to a model translation 
of the J 8 i 4 -band counts (Gardner 1998). Figure 4 shows number counts segregated by 
morphological type. For comparison, the number counts by type from Abraham et 
ah (1996b) and Glazebrook et ah (1995) are shifted to the Hi q hlter and plotted. There is a 
clear difference in relative depth between the optical and IR data for different morphological 
types. For the redder ellipticals, the NIG2 counts reach as deep as the /814 data, but for 
the blue irregulars the data reach only a magnitude shallower than the /814 counts. 
We see good agreement between the Hi q and /814 number counts in each morphological 
bin. In particular, though there are slightly fewer irregular galaxies at Hi q< 20 in the 
6 -band, they are in excellent agreement at fainter magnitudes (the difference is likely 
due to the small numbers of galaxies being counted). Glearly, if we were seeing a strong 
morphological K-correction, this difference would not be the case. Such an effect would 
require that Hi q observations hnd less irregulars at faint magnitudes than predicted based 
on J 814 observations. 

This result implies that the irregular galaxies observed in the MDS are at intermediate 
redshifts {z ~ 0.5). At that redshift, the observed / 814 -band is approximately the rest-frame 
V band, while the observed i/i.e-band is the restframe J-band. If the irregular galaxies were 
at high redshift {z > 2), where optical observations sample the rest-frame ultraviolet, then 
there would be a large morphological difference in the observed near infrared. For example, 
many Lyman Break Galaxies in the HDF are seen to have highly disturbed morphologies 
(see Lowenthal et ah 1996), which have been suggested by O’Gonnell (1997) to be similar 
to the UV morphology of typical spirals. This is consistent with the results of Glazebrook 
et al. (1998), who obtained spectra of 31 galaxies from the MDS irregular galaxy sample 
and found them to be at typical redshifts of 0.5, though they note that the irregulars are 
not necessarily a single uniform population. 

The relatively small number of irregulars in our sample is due to their blue color. 
Glazebrook et al. (1998) hnd irregular galaxies at /814 < 22 have (observed) I-K colors 
of 2-3, with some red outliers. Since the bulk of the irregulars are at I> 24, that would 
imply infrared magnitudes on the order of H>22. This prediction is consistent with our 
data, since we see some irregulars (possibly the red outliers) but a lower percentage than 
was seen in the / 814 -band data. 

Gould we see the high-redshift population in the NIGMOS data? Typical Lyman 
Break Galaxies have IR magnitudes of K=21-22 (Steidel et al. 1996), and may have 
22 < Hi q < 23 and therefore be just at the edge of our ability to classify in the current 
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data {Hi q< 23). Future NICMOS observations should be able to detect such galaxies, 
though much of that work will be done using NIC3 with its coarser plate scale, making 
morphological classihcation more difficult. 

In summary, we have examined NICMOS parallel imaging data. We hud good 
agreement between the Hi q number counts and those in the deepest ground based surveys. 
We have morphologically classihed the galaxies in the NIC2 data following the quantitative 
classihcation scheme of Abraham et al. (1996a). We see evidence for the same population of 
faint blue irregular galaxies observed in the MDS and HDF data. We argue that that the 
irregulars (at i/i.e ~22) are at z< 1. 
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Dehnition Team (STIS IDT) for their encouragement of this project. We thank Jennifer 
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thank P. MacCarthy, L. Yan and L. Storri-Lombardi for useful discussions. We thank the 
referee, Matthew Bershady, for many helpfnl comments which improved the manuscript. 
This work was supported by the STIS IDT through the National Optical Astronomical 
Observatories and by the Goddard Space Flight Center. 


REFERENCES 

Abraham, R.G., Valdes, F. Yee, H.K.C., & van den Bergh, S., 1994, ApJ, 432, 75 

Abraham,R.G., van den Bergh,S., Glazebrook,K., Ellis,R.S., Santiago,B.X., Snrma,P., & 
Griffiths,R.E., 1996a, ApJS, 107, 1 

Abraham, R.G., Tanvir, N.R., Santiago, B.X., Ellis, R.S., Glazebrook, E., & van den Bergh, 
S. 1996b, MNRAS, 279, L47 

Abraham, R.G., 1998, in preparation. 

Bershady, M.A., Lowenthal, J.D., & Koo, D.C., 1998, ApJin press 
Berlin, E. & Arnonts, S., A&AS, 1996, 117, 393 

Bushouse, H., 1997, “The NICMOS CALNICA and CALNICB Pipelines” in The 1997 HST 
Calibration Workshop, eds S. Casertano et al. 

Djorgovski, S. et al. , 1995, ApJL, 438, L13 

Driver, S.P., Windhorst, R.A., Ostrander, E.J., Keel, W.C., Griffiths, R.E., & Ratnatnnga, 
K.U., 1995, ApJ, 449, 23 



Driver, S.P., Fernandez-Soto, A., Couch, W.J., Odewahn, S.C., Windhorst,R.A., 
Philhpps,S., Lanzetta,K., & Yahil, A., 1998, ApJL, 496, L93 

Ellis, R.S., 1997, ARAA, 35 

Frei, Z., Guhathakurta, P., & Gunn, J.E., 1996, AJ, 111, 174 
Gardner, J.P., 1998, PASP, 110, 291 

Gardner, J.P., Cowie, L.L., Wainscoat, R.J., 1993, ApJ415, 9 
Glazebrook, K., Ellis, R., Santiago, B., & Griffiths, R., 1995, MNRAS275, L19 
Glazebrook, K, Abraham, R., Santiago, B., Ellis, R., & Griffiths, R., 1998, MNRAS in press 
Griffiths, R. E., et ah 1994, ApJ, 435, L19 

Holfeltz, S.T., MacKenty, J.W., Sparks, W.B., & Axon, D., 1998, BAAS 191.3.09 
Howell, S.B., 1989, PASP, 101, 616 
Hutchings, J.B., 1995, AJ109, 928 
Kron, R. G. 1980, ApJS, 43, 305 

Loveday, J., Peterson, B.A., Efstathiou, G., Maddox, S.J., 1992, ApJ, 390, 338 
Lowenthal, J.D. et ah , 1997, ApJ, 481, 673 

Moustakas, L.A., Davis, M., Graham, J.R., Silk, J., Peterson, B.A., & Yoshi, Y., 1997, ApJ, 
475, 445 

O’Connell, R. W., 1997, in The Ultraviolet Universe at Low and High Redshift, p.ll, eds. 
Waller, W.H. 

Oke, J.B. 1974, ApJS, 27, 21 

Stanford, S.A., Eisenhardt, P.R.M., & Dickinson, M., 1995, ApJ, 450, 512 

Steidel, C.C., Giavalisco, M., Pettini, M., Dickinson, M., & Adelberger, K.L., 1996, ApJ 
Letters 462, L17 

Thompson, D., Djorgovski, S, & Trauger, J., 1995, AJ, 110, 963 

Thompson, R. L, Rieke, M., Schneider, G., Hines, D. C., & Corbin, M., R. 1998, ApJ, 492, 
L95 

Tyson, A., 1988, AJ96, 1 
Williams et al 1996, AJ, 112, 1335 

Yan, L., McCarthy, P. J., Storrie-Lombardi, L. J., & Weymann, R. J. 1998, ApJ, in press 
Yoshii, Y. & Takahara, F., 1988, ApJ, 326, 1 



- 9 - 


Fig. 1.— Number counts for the Hi q filter. The blue circles show NIC2 and the darker 
blue sqaures show NIC3. Arrows indicate la upper limits. The red + symbols indicated 
the ground-based K-band data (as compiled by Gardner, 1998 and Bershady et ah, 1998) 
translated to the NICMOS Hi q hlter. 

Fig. 2.— Morphological classihcations of galaxies in the NIC2 fields. Classifications by eye 
are indicated by the symbols (stylized asterisks for irregnlars, ellipses for E/SO, and stylized 
spirals for later than SO). The data are plotted in the plane of central concentration (C) vs. 
asymmetry (A). The lines indicate the division of morphological types according to Abraham 
et ah (1996). 

Fig. 3.— A montage of classifiable galaxies from the NIC2 images. From top to bottom, 
galaxies are arranged in bins of ellipticals, spirals and irregnlars according to the antomated 
classihcation. Within each bin galaxies are arranged in descending order of apparent 
magnitnde. Each image is 3.1"x 3.1". Objects identihed differently by eye (see hgnre two) 
are indicated by the letter E,S, or I. 

Fig. 4.— Nnmber Connts by morphological type. The hlled circles are the nnmber connts 
for galaxies that can be morphologically classified. The open sqnares indicate nnmber connts 
of Abraham et ah (1996b) and Glazebrook et ah (1995) shifted to the Hi ^ hlter. 
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